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ABSTRACT 

A physically consistent model of magnetic field generation by convection in a rotating spherical shell with 
a minimum of parameters is applied to the Sun. Despite its unrealistic features the model exhibits a number 
of properties resembling those observed on the Sun. The model suggests that the large scale solar dynamo is 
dominated by a non-axisymmetric m = 1 component of the magnetic field. 
Subject headings: Convection - Magnetohydrodynamics - Sun: dynamo - sunspots 



1. INTRODUCTION 

During the past two decades the understanding of the pro- 
cess of magnetic field generation by motions of an electrically 
conducting fluid has progressed rapidly and increasingly real- 
istic models of dynamos operating in planetary interiors have 
been achieved. Models of the geodynamo, of the ancient Mar- 
tian dynamo or of Mercury's dynamo are primarily restricted 
by the lack of knowledge of certain properties of the respec- 
tive planetary core. The dynamo process operating in the Sun, 
on the other hand, is still subject to controversies with several 
competing proposals. The density variations of several orders 
of magnitude in the solar convection zone and the complica- 
tions introduced by the compressibility of the fluid have pre- 
vented so far a fully convincing model of the solar cycle. 

In this situation it can be helpful to consider the properties 
of a most simple physically consistent model of convection 
driven dynamos in a rotating spherical shell in order to see 
which properties might resemble those of the solar cycle. As- 
suming a given radius ratio of the spherical fluid shell the min- 
imum number of dimensionless parameters is four, namely 
the Rayleigh number R, the rotation parameter r, the Prandtl 
number P denoting the ratio of viscous and thermal diffusiv- 
ity, and the magnetic Prandtl number P„, denoting the ratio of 
viscous and magnetic diffusivity. 

The success of numerical simulations of the generation of 
magnetic fields by turbulent fluid motions relies of the con- 
cept of eddy diffusivities which represent the diffusive actions 
of small scale turbulence that remains unresolved in the nu- 
merical schemes. The ratios of those eddy diffusivities are 
expressed by the values of P and P m used in the numerical 
computations. Frequently the assumption P = P m = 1 is made 
based on the argument that the turbulent diffusion affects the 
transport of heat, momentum and magnetic flux in the same 
way. This assumption is not correct, however, since the trans- 
ports of scalar and vector quantities are aff ected in different 
ways dKays et al.l 12005b lYousef et al.1 120031) . Moreover, the 
diffusivities in the absence of any motion often differ by many 
orders of magnitude, - in particular when the radiative trans- 
port is taken into account in the heat diffusivity -, such that 
even highly turbulent fluids do not fully erase differences in 
the effective diffusivities. 

Even if it is assumed that ratios of the effective diffu- 
sivities for heat, momentum and magnetic flux do not dif- 



fer much from unity in the highly turbulent solar convection 
zone, it must be realized that numerical simulations of con- 
vection driven dynamos depend sensitively on the parame- 
ters P and P,„ near their values of unity. In particular it has 
been found that a transition between two different types of 
dyn amos occurs which is h ysteretic as a function of P and 
P m (ISimitev & Busse 1 120091) . It thus appears to be of inter- 
est to study simple convection driven dynamos in dependence 
on these two parameters in addition to the Rayleigh number 
and the rotation parameter. More realistic solar dynamo mod- 
els involve many more parameters and because of the high 
demand for numerical resolution the opportunities for com- 
prehensive studies of parameter dependences are restricted. 

For the present study a numerical model for convection 
driven dynamos that previously has been employed for inves- 
tigations motivated by problems of the geodynamo (see, for 
example. ISimitev &~B usse 2005) has been modified for solar 
applications. A much thinner spherical shell is assumed and 
no-slip conditions at the inner boundary and stress-free condi- 
tions at the outer boundary are used. The Boussinesq approxi- 
mation has been retained even though it is not appropriate for 
any realistic solar dynamo model and causes an unrealistic 
dependence of the solar differential rotation on depth. Never- 
theless, as we hope to demonstrate, there are properties of the 
dynamo solutions that resemble those observed on the Sun. 

2. MATHEMATICAL FORMULATION 

We consider a spherical fluid shell rotating about a fixed 
vertical axis. The basic equations permit a static solution with 
the temperature distribution 

T s = T + AT V r- l (l-r 1 r\ 

where r denotes the distance from the center of the spherical 
shell measured in multiples of the shell thickness d. The ratio 
of inner to outer radius of the shell is denoted by r). AT is thus 
the temperature difference between the two boundaries. The 
gravity field is given by g = —djr. In addition to d, the time 
d 2 1 v, the temperature v 2 / 7 cm/ 4 , and the magnetic flux density 
v(\i.Q) x l 2 jd are used as scales for the dimensionless descrip- 
tion of the problem where v denotes the kinematic viscosity 
of the fluid, k its thermal diffusivity, q its density, and /i its 
magnetic permeability. The equations of motion for the ve- 
locity vector u, the heat equation for the deviation from the 
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Figure 1, (Color online) Convection-driven dynamos as a function of R, P and P„, for r = 2000. Decaying dynamos are indicated by black dots, MD dynamos are 
indicated by blue diamonds, FD dynamos are indicated by red stars, quadrupolar dynamos are indicated by green triangles, and coexisting MD and FD dynamos 
are indicated by pink squares. 



static temperature distribution, and the equation of induction 
for the magnetic flux density B are thus given by 

<9,u + u- Vu + rkxu = -V7r + 6r + V 2 u + B- VB, (1) 
Vu = 0, (2) 

P(d t e + u ■ V9) = (Rr]r-\l-r)T 2 )r ■ u + V 2 9, (3) 
V-B = 0, (4) 

V 2 B = P m (d,B + u-VB-B-Vu), (5) 

where d t denotes the partial derivative with respect to time t 
and where all terms in the equation of motion that can be writ- 
ten as gradients have been combined into Vtt. The Boussi- 
nesq approximation is assumed in that the density g is re- 
garded as constant except in the gravity term where its temper- 
ature dependence, given by a = -(d g/ dT) / g =const., is taken 
into account. The Rayleigh numbers R, the Coriolis number 
r, the Prandtl number P and the magnetic Prandtl number P,„ 
are defined by 



a-fATd 4 



2fld 2 V V 
p = — p = — 

V K A 



(6) 



where A is the magnetic diffusivity. Because the velocity field 
u as well as the magnetic flux density B are solenoidal vec- 
tor fields, the general representation in terms of poloidal and 
toroidal components can be used 



u = V x (Vv x r) + Vw x r 
B = Vx(V/ixr) + Vgxr 



(7) 
(8) 



Equations for v and w are obtained by multiplication of the 
curl 2 and of the curl of equation (1) by r. Analogously equa- 
tions for h and g are obtained through the multiplication of 
equation (5) and of its curl by r. 

No-slip boundary conditions will be used at the inner 
boundary and stress-free conditions are applied at the outer 
boundary, while the temperature is assumed to be fixed at both 



boundaries, 



v = d r v = w = = at r = r; = fy/( 1 — rj) 



(9) 



v = «9 2 v=<9,.(w/r) = e = atr = r = 1/(1-77). (10) 

For the magnetic field, electrically insulating boundaries are 
assumed such that the poloidal function h must be matched to 
the function which describes the potential fields outside 
the fluid shell 

g = h-h (e) = d r (h-h (e) ) = at r = n and at r = r . (11) 

Throughout this paper we shall use i] = 0.65. 

The numerical integration proceeds with a pseudo-spectral 
method as described by Tilgner (1999) which is based on an 
expansion of all dependent variables in spherical harmonics 
for the 9, ^-dependences, i.e. 

v = ^V, m (r,/)P, m (cos0) exp(im<p) (12) 

and analogous expressions for w, Q,h and g. Here P'" denotes 
the associated Legendre functions. For the r-dependence ex- 
pansions in Chebychev polynomials are used. For the compu- 
tations to be reported in the following a minimum of 41 col- 
location points in the radial direction and spherical harmonics 
up to the order 128 have been used. 

The numerical mod el outlined above i s essen tially the same 
as has been used by iGilm an & M iller I dl981l) and iGilman I 
d 1 983b . At that time, 30 years ago, computer resources al- 
lowed time integrations only for a few selected parameter val- 
ues and the phenomenon of bistability (see below) was not 
yet known. It thus seems appropriate to return to the most 
simple physically consistent formulation of convection driven 
dynamos in rotating spherical fluid shells as formulated in this 
section. 

3. DYNAMO SOLUTIONS IN DEPENDENCE ON P AND P m 

In choosing parameter values in our model we follow 
IGilman I ( 11983b in selecting a value around 2000 for r which 
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Figure 2. (Color online) Coexisting nonlinear dynamo solutions at identical 
parameter values in the case r) = 0.65 r = 2000, R = 150000 and P/P,„ = 0.2. 
The fluctuating poloidal magnetic energy component M p (solid circles) and 
the mean poloidal magnetic energy component M p (empty circles) are scaled 
on the left ordinate, while the Nusselt number at r = r, (solid diamonds) is 
scaled on the right ordinate. Cases started from initial conditions on the FD 
branch are shown in red and connected by dashed lines, and those started 
from initial conditions on the MD branch are shown in blue and connected 
by dotted lines. 



in the case of the Sun corresponds to an eddy viscosity v = 
1 .6 x 10 8 m 2 /s. Values for R have usually been selected within 
a range of up three times the value needed for dynamo action. 
Otherwise the dynamos become too chaotic and regular fea- 
tures disappear. Our study has thus focussed on the variations 
of dynamos with the parameters P and P m . 

In figure 1 the existence of solutions for convection driven 
dynamos in the (P,P,„,/?)-space is outlined. In order to 
achieve dynamo action the Rayleigh number R must be suffi- 
ciently high such that the magnetic Reynolds number R m ex- 
ceeds a critical value. The latter is inversely proportional to 
the magnetic Prandtl number P,„. Since our study focusses 
on reasonably regular dynamos we have found it difficult to 
obtain dynamos for small values of P,„ since dynamos will 
become highlychaotic in this case or even decay. As indi- 
cated in figure[T]there exist a large variety of dynamos ranging 
from dipolar to quadrupolar ones. At higher Rayleigh num- 
bers the dynamos are often of mixed polarity. While most of 
the dynamo solutions have been obtained for the radius ratio 
r) = 0.65, some cases obtained for r\ = 0.6 have been added in 
the figure since it was found that the dynamos depend only 
weakly on this parameter. 

Of particular interest is the distinction between MD- 
dynamos which are dominated by a strong and nearly steady 
axisymmetric dipolar magnetic field and the FD-dynamos for 
which magnetic energies of the axisymmetric components are 
far less than those of the non-axisymmetric components and 
which exhibit a cyclic behavior. Over a considerable extent 
of the parameter range these two types of dynamos exist si- 
multaneously at the same values of the external parameters 
as demonstrated by the example of figure [2] It is remarkable 
that the two types of dynamos exhibit nearly the same Nusselt 
number. This property indicates that both types are similarly 
efficient in the use of the available buoyancy. For more de- 
tai ls on the phenomenon of bistability we refer to the paper 
bv ISimitev & Busse I d2009l) . A comparison of the magnetic 
and kinetic energy densities for the two types of dynamos is 
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Figure 3. (Color online) Coexisting distinct nonlinear dynamo solutions at 
identical parameter values - a MD (left column (a,b)) and a FD dynamo 
(right column (c,d)) both at R = 1.5 X 10 6 , r = 2 X 10 4 , P = 1 and P,„ = 5. 
Panels (a,c) and (b,d) show time series of magnetic dipolar energy densities 
and kinetic energy densities, respectively. The component X p is shown by 
thick solid black line, while X t , X p , and X t are shown by red, green and blue 
lines, respectively, and labeled. X stands for either M or E. 

shown in figure [3] Here the definitions 

M p = (| V x (V/z x r) | 2 }/2, M, = <| Vg x r | 2 )/2 (13) 

M P = (| V x (V/z x r) | 2 }/2, M t = (\ Vg x r | 2 )/2 (14) 

have been used where the angular brackets indicate the aver- 
age over the fluid shell and where h refers to the azimuthally 
averaged component of h and h is given by h = h - h. Anal- 
ogous relationships hold for the kinetic energy densities for 
which M,h,g are replaced by E, v, w. 

Near the end points of the bistability regions transitions 
from MD- to FD-dynamos and vice versa may, of course, 
occur owing to exceptionally large fluctuations. Such a tran- 
sition was observed in the long-time dynamo integration of 
iGoudard & Dormvl J2008I) . As their convection driven dy- 
namo changed from the MD-state to the FD-state these au- 
thors noted a solar like behavior in butterfly-like diagrams in 
which they had plotted the axisymmetric (^-component of the 
magnetic field, B v , as a function of latitude and time. As must 
be expec ted from the decr ease with depth of their differential 
rotation (I Yoshimural 1 9751) . the waves of the butterfly diagram 
propagated to higher latitudes instead of lower ones. It turns 
out, however, that non-axisymmetric components of the az- 
imuthal magnetic field are comparable in amplitude to the ax- 
isymmetric component B v . It thus appears to be appropriate 
to study these dynamos in more detail. 

In figure |4] a typical example of a nearly periodic FD- 
dynamo is shown. In the four columns of plots the devel- 
opment of various components of the magnetic field is shown 
in five equidistant steps in time covering about half a period 
of the magnetic cycle. As is evident from the radial compo- 
nent B r displayed in the third column, the dynamo appears 
to operate mainly in one of two meridional hemispheres. In 
the second column maxima of the horizontal magnetic field 
at the depth of 0.1 from the surface have been indicated as 
a measure of the likelihood of the appearance of sunspots. 



4 



Simitev and Busse 




Figure 4. (Color online) Approximately half a period of a dynamo oscillation in the case r = 2000, R = 120000, P = 1.2, P m = 4.5. The first column shows 
meridional lines of constant B v on the left and poloidal fieldlines, rsmOdh/ 89 on the right. The second column shows lines of constant horizontal magnetic 
field (B 2 V +B^)'/2 at r = 0.9 corresponding to -0.9, -0.8, -0.7, 0.7, 0.8, 0.9 of its maximum absolute value with the sign determined by the sign of B v , and the 
third column shows lines of constant B, at r = r . The last column shows Re(dg"' =1 / 88) to the left and Im(9g"' =l /86) to the right. The five rows are separated 
equidistantly in time by At = 0.0224. 
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Figure 5. "Butterfly" diagrams: B™ =0 + |B™ =1 |sgn(fi™ =0 ) (top) and B" r ' =0 + 

|2?"' =1 |sgn(S" !=0 ) (bottom) are plotted as functions of latitude and time in the 
case P = 1.2, r = 2000, R = 120000, P„ = 4.5. 

Since the underlying hemispherical structure hardly moves 
the pattern of the second column seems to reflect the phe- 
no menon of the active lo ngitude of sunspots (see, for exam- 
ple, lUsoskin et al.l d2007l) and r eferences therein) or of solar 
X-rav flares dZhang et al.ll2007l) . Since sometimes the m = 2- 
component of the magnetic field dominates instead of the 
m = 1 -component the phenomenon of two active longitudes 
180 degrees apart can also be observed in the dynamo sim- 
ulations. Persistent active longitud es have been detecte d in 
various types of active stars as well (lUsoskin et al.ll2007h . 

In figure[3]two "butterfly" diagrams have been plotted. Be- 
cause the strength of the axisymmetric component of B v alone 
is not representative of the horizontal magnetic field that may 
give rise to sunspots we have added in the upper plot the main 
non-axisymmetric component characterized by the azimuthal 
wavenumber m = 1 . While there is some wave-like movement 
towards lower latitudes, the main pattern moves towards high 
latitudes. This must b e expected s i nce th e differential rotation 
decreases with depth lYoshimural ( 11975b . The movement to- 
ward higher latitudes is even more pronounced in the pattern 
of the radial component B r of the magnetic field as shown in 



the lower plot of figure [5] A phase shift between the patterns 
of the upper and the lower plot can also be noticed in that the 
sign change of B r occurs after that of B v . This feature agrees 
with the well known solar phenomenon that the change in the 
sign of B r occurs near the maximum of the sunspot cycle. 

4. CONCLUSION 

Most models of the solar cycle are based on essentially ax- 
isymmetric dynam os, but there are exceptio ns such as the dy- 
namo proposed by iRuzmaikin et alJ (119881) based on a mean 
field approach. It is not quite clear how strong the non- 
axisymmetric components of the solar toroidal magnetic field 
with wavenumbers m = 1 or m = 2 are in relation to the ax- 
isymmetric comp onent. Estimates ra nge from 0.1 based on 
sunspot numbers (lUsoskin et a l. 2007) to 0.43 based on pow- 
erful X-ray flares (IZhang et al.ll2007l) . A main conclusion to 
be derived from the dynamo model outlined in this paper is 
that the large scale non-axisymmetric components are likely 
to be an important ingredient of solar magnetism. The well 
documented existence of active longitudes certainly supports 
this conclusion. 

There are, of course, also major differences between so- 
lar phenomena and certain features of the present Boussinesq 
model. Some of these are well understood as for example the 
difference in the propagation of dynamo waves. This is most 
probably caused by the fact that the model always exhibits 
a differential rotation decreasing with depth for FD-type dy- 
namos, while the Sun possesses a differential rotation that in- 
creases with depth in the uppermost 15% of the convection 
zone. For a preliminary attempt to take this effect into account 
we ref er to a recent companion paper by ISimitev & Busse I 
JIOH) . 
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